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Abstract 



We discuss a hydrodynamic obstruction to bubble wall acceleration during a cos- 
mological first-order phase transition. The obstruction results from the heating of the 
plasma in the compression wave in front of the phase transition boundary. We provide 
a simple criterion for the occurrence of the obstruction at subsonic bubble wall velocity 
in terms of the critical temperature, the phase transition temperature, and the latent 
heat of the model under consideration. The criterion serves as a sufficient condition of 
subsonic bubble wall velocities as required by electroweak baryogenesis. 



1 Introduction 



Cosmological first-order phase transitions can lead to many interesting phenomena during 
the evolution of the early universe, such as electroweak baryogenesis [I] or the production of 
a stochastic background of gravitational waves [2HI] . A first-order phase transition proceeds 
by bubble nucleation and subsequent bubble expansion and one essential quantity in the 
description of phenomena linked to the phase transition is the velocity of the expanding 
bubble walls £ w . For example, the electroweak baryogenesis mechanism is based on the 
diffusion of particle asymmetries into the plasma in front of the bubble wall and subsonic 
bubble walls are necessary to build up a baryon asymmetry. On the other hand, fast moving 
walls are essential for the production of a sizable amount of gravitational radiation by bubble 
collisions [3H7], turbulence [8j[9] or magnetic fields [9]. 

The analysis of the bubble wall velocity generally assume that after a short period of 
acceleration of order ~ i/T, (where T represents the typical energy scale associated with 
the temperature or latent heat of the transition) the pressure difference that drives the bubble 
expansion is balanced by friction and the bubbles subsequently expand with a constant speed. 
To quantify this friction requires solving a coupled system of Boltzmann equations for all 
particle species with a sizable coupling to the Higgs field. This intricate calculation has so far 
only been performed in the Standard Model (SM) [10] and in the Minimal Supersymmetric 
Standard Model (MSSM) [TT] under the assumption of small wall velocities. 

On the other hand, in the limit of highly relativistic wall velocities it is found that the 
friction in the plasma tends to a constant [12] (up to possible log(7„,) corrections), opening the 
possibility of continuously accelerating (runaway) bubble walls when the pressure difference 
along the phase boundary overcomes this threshold. This runaway behavior is realistic in 
many models, under the assumption that no hydrodynamic obstruction prohibits that the 
highly relativistic regime is reached. 

The goal of this paper is to analyze one of these possible obstructions based on the 
heating of the plasma in front of the phase boundary during bubble expansion. Implicitly, 
this effect was already observed in refs. [10l[T3l[Tl] where finite wall velocities have been 
reported in the limit of vanishing or at least very small friction. In the recent work [T4] this 
result was obtained under the assumption that the temperature in the Higgs wall is identified 
with the temperature in front of the phase boundary. Furthermore, the analysis focused on 
models with an equation of state similar to the standard model. Here, we relax those two 
assumptions and present a simple criterion for the occurrence of the obstruction. If our 
criterion holds in a specific model, the wall velocity is subsonic and electroweak baryogenesis 
is in principle possible. The heating effect only provides an upper limit and a concrete 
determination of the wall velocity still requires some knowledge of friction. However, as long 
as friction is not too strong (as it is e.g. the case in the SM) the resulting wall velocities 
are fairly accurate. Besides, the baryon asymmetry in electroweak baryogenesis is not very 
sensitive to the wall velocity as long as it is significantly below the speed of sound and large 
enough to avoid a saturation of the sphaleron process (see e.g. JT5] or [IS])- Hence, the 
knowledge of the precise wall velocity is not so relevant for baryogenesis as long as it is 
subsonic. 

The paper is organized as follows: in the next section the hydrodynamic treatment of the 
system is reviewed, and we briefly sketch the origin of the obstruction. In sec. |3]the heating 
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effect is discussed analytically in a system with a bag equation of state. We then study how 
the heating effect impacts the acceleration of the wall in a toy model and present a simple 
criterion for the occurrence of the obstruction in sec. HI In sec. [5] we briefly comment on the 
interplay between the hydrodynamic obstruction and friction before we apply our results to 
specific models and conclude in sec. El 



2 Origin of the obstruction 

In this section we introduce the basic concepts and set up the notation used for the hydro- 
dynamic analysis of the system of expanding Higgs bubbles [TT14T9] . For most parts we use 
the notation of [13] . 

The dynamics of the combined "wall-plasma" system is determined by the equations of 
motion of the plasma and of the Higgs field. However, in the present case it is advantageous 
to replace the equation of motion of the plasma (which would be of Boltzmann type) by the 
assumption of local thermal equilibrium and energy-momentum conservation, leading to the 
hydrodynamic approximation. The energy-momentum tensor of the Higgs field <fi is given by 



ld P 4>d^-v ((l>) 



X 



where Vo(0) is the renormalized vacuum potential. If the plasma is locally in equilibrium its 
energy-momentum tensor can be parametrized as 

T^ sma = wu^-g^p, (2) 

where w and p are the plasma enthalpy and pressure, respectively. The quantity is the 
four-velocity field of the plasma, related to the three-velocity v by 

^ = ^== = (7,7v). (3) 

A constant (ft background contributes to the total pressure [see eq. (1)] and from now on we 
will use p for this total pressure, including such contribution. The enthalpy w, the entropy 
density a and the energy density e are defined by 

dp _ dp _ rr dp . s 

W = T dT' a =df' e = T dT- p ' (4) 

where T is the temperature of the plasma. One then has 

w = e + p . (5) 

Conservation of energy-momentum is given by 

dT^ = «9%t + d»Tf™™ = . (6) 

We are interested in a system where the bubble expands at a constant speed and, assuming 
there is no time-dependence, eq. ([6]) reads in the wall frame (with the wall and fluid velocities 
aligned in the z direction) 

d z T zz = d z T z0 = . (7) 
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Integrating these equations across the phase boundary and denoting the phases in front and 
behind the wall by subscripts + (symmetric phase) and — (broken phase) one obtains the 
matching equations (in the wall frame): 

o o 99 9 9 / \ 

w + v + r ) + + p + = W7_t>_7_ +p_ , w+v + y+ = w_w_7_ . (8) 
From these equations we can obtain the relations 



P+-P- v + g- +P+ / n 

V+V- = , — = ; • (9) 

e + — e_ v - e + + p- 

In a concrete model the thermodynamic potentials can be calculated in the two phases 
and the temperature at which the phase transition happens can be determined using the 
standard techniques [21] ■ Still, there are three unknown quantities (T_, v + and t>_) and 
only two equations fl5J), so that up to this point all hydrodynamically viable solutions are 
parametrized by one parameter. We will parametrize the solution by its wall velocity £ w . 
Ultimately, the wall velocity is deduced from the equation of motion of the Higgs 

where fa denotes the particle distribution function of the ith species. By decomposing 

fi(p) = f?(p) + Sfi(p) , (11) 



where f^ q = l/[exp (Ei/T) =F 1] is the equilibrium distribution function with Ef = p 2 + m 
eq. ( TTOT) takes the simple form 



87 

— -/C(0) = O. (12) 

The second term contains the free energy J 7 that drives the expansion of the bubble and 
K,{4>) stands for the friction term that arises from deviations of the particle distributions in 
the plasma from equilibrium. In principle, calculation of /C(</>) involves solving a coupled 
system of Boltzmann equations for all particle species with a sizable coupling to the Higgs 
field. This intricate calculation has been performed in the Standard Model [10] and in the 
MSSM [11] and under the assumption that the deviation from thermal equilibrium is small, 
i.e. 5fi(p) <C fi(p), which is only true for weakly first-order phase transitions. 

The aim of the present work is to show that the bubble wall velocity can approach 
a subsonic value even if the friction is very small - contrary to the naive expectation of 
supersonic velocities or even runaway behavior in the limit of very small friction. The 
occurrence of this effect can be understood by inspecting eq. (|T2|) . Assuming that in the 
steady state the bubble is large enough so that one can use the planar limit, one obtains by 
integration the pressure (in the wall frame) that drives the expansion 



dr 



dJ 7 

dzd z (j) — — = I dzd z 4>K, = Ff r . (13) 



The physical interpretation of this equation is that the change of pressure in the wall drives 
the expansion of the bubble and this driving force F^ r is ultimately balanced by the friction 
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force Ff r in order to reach a stationary state. Without the influence of the bubble, this change 
in pressure is always positive, since nucleation requires the temperature in the symmetric 
phase to be below the critical one. However, some particles are reflected at the phase 
boundary, leading to a heating effect in front of the wall. Hence the temperature experienced 
by the phase boundary is increased and a hydrodynamic obstruction can occur when the 
bubbles accelerate while building up a compression wave in front of the Higgs wall. At some 
velocity, the average temperature^ in the wall might approach the critical one, the driving 
force goes to zero and the bubble cannot further accelerate even in the limit of vanishing 
friction as defined by (fl2]l and ( fl3l) . 

As a final remark in this section, notice that there is a certain freedom in splitting the 
last term of ( JTUi) into equilibrium and friction contributions. At the critical temperature T c 
(where the pressure in the two phases is equal), the only solution to ( 112)) is a static wall. The 
plasma is at rest and the temperatures on both sides of the Higgs wall coincide: the system 
is in equilibrium everywhere. For a temperature below the critical one, all deviations from 
equilibrium in the particle distributions are proportional to the wall velocity and also the 
friction term vanishes in the limit of zero velocity. This leaves a certain arbitrariness in the 
temperature profile of the equilibrium distributions f^ q in (11) along the wall and different 
choices for the temperature profile in the equilibrium distributions will result in different 
friction terms. 

A necessary feature of the temperature profile is that it interpolates between the two 
values T + and TL. This ensures that the deviations from equilibrium 5fi vanish away from 
the bubble wall. Besides, an upper bound on the wall velocity will be derived in the following 
by neglecting the friction term. Hence, one has to ensure that the friction term is always 
positive and reduces the wall velocity for the chosen temperature profile. A very convenient 
choice for the temperature profile is obtained by enforcing conservation of the two energy- 
momentum tensors resulting from the particle distributions f^ q and 5fi separately. First, 
this has the advantage that in the limit of infinite interaction rates the particle distribution 
functions f^ q are the physical ones and 5fi will vanish. Second, for this choice entropy 
is conserved when the deviations 5fi are neglected. A simple phenomenological approach 
for the friction term shows (see e.g. ref. [22]) that negative friction terms would lead to 
entropy decrease. This indicates that in general the H-theorem ensures that the expansion 
around this equilibrium leads to positive friction terms that reduce the wall velocity. We 
will essentially utilize this choice for the temperature profile in section HJ 



In this section we discuss a system whose equation of state is in the broken phase given by 



where e denotes the false-vacuum energy resulting from the Higgs potential, while in the 
symmetric phase 



3 The heating effect 



v- = 




aJT_ — e , 



(14) 



P+ = 3 a + T + > e + = a + T + > 



(15) 



We will be more specific what the average temperature is in sec. 2] 
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with a different number of light degrees of freedom across the wall and therefore different 
values a + and a_ (with a + > a J) and different temperatures on both sides of the wall. 
These expressions correspond to the so-called bag equation of state. This approximation 
works reasonably well when the Higgs vev does not change much in between the critical 
temperature and zero temperature and particles can be independently of the temperature 
divided into "light" and "heavy" in the two phases. In particular, this is the case for models 
where the potential barrier between the broken and symmetric phases is present even at 
vanishing temperature, generally leading to relatively strong phase transitions. 

For temperatures close to the critical one, it is more physical to parametrize the pressure 
and energy density in the symmetric and broken phases by 

P-^\a + Ti-t c (T-/T c -l) , P+ ~~a + T±, (16) 

e_~a+Ti-e c , e + ~a + T^, (17) 

and £ c is the latent heat at the critical temperature. Comparing with the bag equation of 
state one can identify 

4 = ^(a+-a_)T c 4 = 4e. (18) 

Our analysis will be done mostly assuming the bag equation of state in order to facilitate 
comparison with the existing literature. However, when rephrased in terms of the latent 
heat all results are equally applicable to models with a weak phase transition where the bag 
equation of state does not hold. 

Using the bag equations of state ( fi4|) and ( 1T51) in eq. (jHJ) we get 



1 - (l-3a. 



3-3(l + a+)r 
v + 3 + (1 — 3a + )r 
v- l + 3(l + a+)r 



(19) 



where we defined 



a + = — -j = - — —j , r = — ^ . (20) 

The quantity a + is the ratio of the vacuum energy to the radiation energy density and 
typically characterizes the "strength" of the phase transition: the larger a + the stronger the 
phase transition. These two equations can be combined to give 



1 




(21) 



so that there are two branches of solutions, corresponding to the ± signs in eq. ([211 . 

Generally, there are three different types of solutions for the bulk fluid motion (recent 
reviews on the topic can be found in [14J and [20]): Detonations, deflagrations and hybrid 
solutions. In detonations, the wall expands at supersonic velocities and the vacuum energy 
of the Higgs leads to a rarefaction wave behind the bubble wall, while the plasma in front of 
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the wall is at rest. In this case, the wall velocity is £ w = v + > u_, and therefore detonations 
are identified with the + branch of solutions in (l2~Tj) . In deflagrations, the plasma is mostly 
affected by reflection of particles at the bubble wall and a compression wave builds up 
in front of the wall while the plasma behind the wall is at rest. In this case, the wall 
velocity is identified with £ w — t>_ > v+, corresponding to the — branch of solutions in ( I2TT) . 
"Pure" deflagrations are subsonic, while the hybrid case occurs for supersonic deflagrations 
where both effects (compression and rarefaction wave) are present. In the following we are 
concerned with deflagrations. 

In the following we present a rough constraint for the occurrence of the obstruction. We 
assume that the bag equation of state holds even locally in the wall such that the pressure 
is given by the simple expression p = |a(0)T 4 . Then, the driving force ( ITS]) reads 

dzd z <f> ^- = e- \ [ daT 4 . (22) 



dr 
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For a monotonously changing temperature in the wall, the last contribution is bounded by 
(we will see later on that T_ < T+) 

^(a + - a_)Tt < e - F dr < i(o+ - a_)Tl . (23) 

Analogously to a + one can define the quantities a^v, (*- and a c as the vacuum energy 
normalized to the energy density in the symmetric phase, behind the wall and at the critical 
temperature 

a - = Z5v aN = ^h- (24) 

«„= « =l5tZ5z<l (25 ) 
a + T c 4 3 a + 3 

The nucleation of bubbles happens at temperatures below the critical one, a c < a^. After 
nucleation, the bubble wall accelerates, increasing T + and therefore decreasing a+. Likewise, 
T_ is expected to increase and a_ to decrease during the acceleration. If the phase transition 
is sufficiently weak, this heating effect eventually leads to a vanishing driving force and the 
acceleration of the wall ceases. According to (123!) this has to happen somewhere in the 
window 

T_ < T c < T + . (26) 

We will show that for small wall velocities the relation T + > T_ > Tjy holds, so the hydro- 
dynamic obstruction definitely occurs in models where T/v is very close to T c . 

In the following we briefly review the solutions of the plasma velocity in the limit of 
small wall velocities £ w as done in [13]. For deflagrations, the fluid in the bubble is at rest 
and changes in the bubble wall to a finite value according to (j2T|) where v _ is identified 
with the wall velocity £ w . The fluid motion in the compression wave is then determined by 
hydrodynamics. Since there is no intrinsic macroscopic length scale present in the system, 
solutions to the hydrodynamic equations are self-similar and only depend on the combination 
£ = r/t, where r denotes the radial coordinate of the bubble and t the time since nucleation. 
The plasma then fulfills the equations 

' dp, (27) 



2| =7 2 (l-<) 



^-1 

n2 
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and 



w 

where c s = l/v3 denotes the velocity of sound in the plasma and /i(£,t>) is the Lorentz- 
transformed fluid velocity 

K^v) = ^-. (29) 
An example of the fluid motion and the enthalpy in case of a deflagration is shown in Fig. [TJ 
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Figure 1: Example for a deflagration. The left plot shows the fluid velocity v(£) while the right 
plot shows the enthalpy w(£). 



The boundary conditions for deflagrations read in the rest frame 

v w = v(£ w ) = fi(v-,v+) , 



(30) 



(31) 



and the fluid velocity in the compression wave is for small wall and fluid velocity t)<C 1] 
given by 

a 

The enthalpy changes in the compression wave between the Higgs wall and the plasma in 
front of the shock by a factor (there is also a jump in the enthalpy in the shock front which 
is however of order 



rp4 
1 N 



+ _ aN _ W + 



cxp 



v) dv 



v(c B ) 



exp (8 £ w v v 



(32) 



So in this limit indeed T + > T N , as required for an obstruction ff26l) . Using the bag equation 
of state and fl2Tj) one finds 

— ~ 1 - 3a + - 6a + £^ , (33) 
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and hence 

v w « 3a+£ w . (34) 

If the obstruction would occur as soon as T + exceeds the critical temperature, the wall 
velocity would be 

a = ^ • (35) 

Nevertheless, T_ is typically smaller than T+, so a more conservative bound arises from the 
pressure in the broken phase, that is determined using 



such that 



and hence (using 



T, a_ w + 1 — 3a 



(36) 



T_ a + W- 1 — 3a 



c - (1 + 120+^) ~ 1 + 12a + ^ , (37) 



log — ~ 3£> + . (38) 

Assuming that the obstruction occurs for T_ = T c one finds analogously to ( |35|) (assuming 
a c ~ a + in a first approximation) 

and hence a slightly larger wall velocity. The obstruction should hence be somewhere in the 
range 

^(3±1)^. (40) 

However, even if this criterion leads to a supersonic wall velocity, the obstruction might 
still be present, since for sizable wall velocities the above approximation underestimates the 
effect of the heating. In order to show this last statement, we assume a large but subsonic 
wall velocity < c s . When the wall velocity becomes supersonic, the expansion mode 
by deflagration is prone to decay into a detonation and we assume that the obstruction is 
absent in this case [T3] . In the following we derive some analytic results under the assumption 
a + < 1 and then present some numerical results for £ w = c s and arbitrary aq_. In this case 
we obtain for the plasma velocity in front of the wall (using (|2Tj) ) 

1 - a/o;+(2 + 3a + ) 2 

v + ~c s ^— p — + 0(^ w -c s f). (41) 

1 + a + 

Notice that positive velocities imply that a + < 1/3. If a + is larger than this value, defla- 
grations are not possible and the bubble expansion proceeds by supersonic detonations. In 
the wall frame the plasma has in front of the wall the velocity 

3 



* - - c 8 ) + y/2a+/3) , (42) 



8 



which is obviously only a valid approximation for (c s — £ w ) < a/2« + /3. Comparing 
with (|4*2l) shows that for wall velocities close to the speed of sound the heating effect in front 
of the wall is much more efficient than in the limit of small wall velocities. The bound (j40p 
is hence a conservative one. 

In the compression wave the plasma velocity decreases to a value v s h and then drops to 
zero in the shock front. The decrease in enthalpy in the compression wave is of order (using 

m) 

log -t ~ Cs (v w - v sh ) + 0{v 2 w , v 2 sh ) , (43) 



T, 



sh 



whereas the jump in enthalpy in the shock front is of order (using (jSj) in the shock front) 



rp4 

-f ~ exp(4c.v. ft ) + 0(v 3 sh ) 

2 N 



We therefore obtain 



log 



N 



c s ) + V2a+/3 



(44) 



(45) 



which is in leading order independent of the value of the fluid velocity at the shock front v s h 
and the fluid profile in the compression wave. 




Figure 2: The heating effect assuming a wall velocity of the speed of sound, £ w = c s . The plot 
shows log T+/T/v as a function a+. The dashed line is the leading order result (|46p while the 
straight line is the fit (|47p to the numerical result. 



To derive a simple criterion for the occurrence of the obstruction, we assume again that 
the average temperature in the wall is close to T + and that supersonic deflagrations can 
decay into detonations. For a fixed value of a + , one can then determine from ( 145)) in leading 
order in a + 

\og^±^y^/2. (£ w = c s ) (46) 

-LN 
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The full numerical result is shown in Fig. [21 and a formidable fit is given by 

log jpt ~ _ A a+ _ i a f = n( a+ ) . = c.) (47) 

If the obstruction would occur as soon as the temperature in front of the wall T + surpasses 
the critical temperature T c , an easy and model independent criterion for the occurrence of 
the obstruction would be 

\og^<n(a c ). (48) 

However, the enthalpy in the bubble interior is in leading order (which is 0{\fa)) given by 

iu_ « (1 - Ac s v w )w + « w N , (49) 

and so the temperature in the bubble is close to the nucleation temperature. In the limit a c — > 
1/3, the temperature behind the wall even vanishes. Hence, one expects that the temperature 
T + has to surpass the critical temperature T c significantly to obtain an average temperature 
in the wall that is close to the critical one. To obtain the wall velocity corresponding to the 
hydrodynamic obstruction (if present) requires the knowledge of the temperature profile in 
the wall in order to determine the average temperature experienced by the Higgs. This is 
the topic of the next section. 

4 Local equations of motion in a toy model 

In the last section, we presented results for the heating effect in front of the wall during the 
bubble expansion. In particular, we exemplified how to calculate the temperature on both 
sides of the Higgs wall, T + and T_, as functions of the wall velocity the strength of the 
phase transition a+ and the nucleation temperature T/v- The topic of the present section is 
to connect this heating effect with the occurrence of the obstruction to bubble acceleration. 
An exact treatment of this question would require to determine the particle distribution 
functions by solving a system of Boltzmann equations. The thermodynamic potentials in 
the wall could then be obtained via (|10p . On the other hand, the obstruction has to happen 
in the window given by eq. (|26|) and any simplifying assumption should lead to reasonable 
bubble wall velocities. In this section, we assume that the ansatz of local equilibrium (121) 
also applies in the bubble wall. This approximation becomes exact in the limit that the 
mean free path of the particles in the plasma is much shorter than the bubble wall thickness. 
Even though the mean free path of some particles is much larger than the wall thickness, we 
expect that the free energy should not depend too much on the out-of-equilibrium features 
of the particle distribution functions since energy-momentum conservation generally holds 
(see also the comments at the end of section [2]). 

4.1 Solving the local equations of motion 

To obtain the steady-state profiles across the phase transition wall of quantities like the 
velocity, temperature and Higgs field, we integrate the following system of coupled differential 
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equations in the planar limit: 




d z [u^ 2 v] = , 
d z ^(^) 2 + w 7 V+p =0 



(51) 



(50) 



(52) 



The first equation is the Higgs equation ([12"]) and to be specific we use the friction term 
(see [H] and also [23] for a similar approach) 



The particular form of the friction term is not relevant since we are interested in the regime 
of small friction, r\ — > 0. The two extra equations correspond to the differential (and static) 
form of energy- momentum conservation They result from the assumption that the 

system is locally in equilibrium, also on scales that are comparable with the wall thickness. 
Their integration across the wall gives immediately Steinhardt's matching conditions ([8]). 

In order to solve the system of equations, we proceed as follows. The basis quantity 
is the free energy of the system J-"(0, T) from which all the thermodynamic potentials can 
be derived. We determine the nucleation temperature of the system using the standard 
techniques of the bounce solution (21]. With a specific wall velocity £ w and this information, 
the hydrodynamic system ([9]) can be solved. Finally, we solve the system of differential 
equations fl50|) - fl52|) and determine which friction coefficient rj corresponds to the wall velocity 
£ w we have chosen. When the wall velocity is increased, the friction coefficient eventually 
turns negative. This signals the occurrence of the hydrodynamic obstruction. The numerical 
results obtained by this procedure depend neither on the bag equation of state nor on the 
assumption of small fluid or wall velocities. 

The system of equations can be solved in the following way: First, equation (l5Tj) can be 
used to obtain the fluid velocity as a function of the Higgs vev and the temperature 



K, = -T N r]vd z (p . 



(53) 




(54) 



Then eq. (I52"j) can be used to obtain a closed expression for d z 4>(4>, T) 



h&M, T)) 2 = w<yV + p- u +1 2 + v 2 + - p + . 



(55) 



Finally, the Higgs equation ( 1501) can be recast using (1521) as 



d z <p 1 ^-d z <piC = d z {uj 1 2 v 2 +p) , 



(56) 



which leads to 



d(f) driui^v 2 + p) 



(57) 



dT K, - d^u^v 2 ) ' 
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For a fixed wall velocity £ w , one unique choice of r/ leads to the correct boundary conditions 
<j){T±) = (j)±. The wall velocity can then be varied to lead to the velocity corresponding to 
the obstruction where rj = 0. 

As a remark we note that solving the system of local equations ( l52l) also allows to 
determine the wall thickness. Usually the wall thickness is determined for the nucleated 
bubbles from the bounce solution. In principle the wall thickness during the bubble expansion 
is expected to be smaller than this value since the acceleration of the wall also compresses 
the wall. In practice we find that this effect changes the wall thickness by only a few percent. 

4.2 A concrete (toy) model 

Consider the effective potential 

V{<f>, T) = A(0 2 - v 2 ) 2 - ?± T 4 + Cl T 2 2 - c 2 T0 3 . (58) 

This model is inspired by systems where the strength of the phase transition comes from 
thermal contributions of bosons coupled to the Higgs that add to the thermal potential terms 
proportional to T<f) 3 . This is for example the case in the light stop scenario [21] of the MSSM 
where the W-bosons and the stops give rise to the cubic terms. Another example of this 
class of models are hidden sector models [25]. The parameters A, c\ and c 2 can be exchanged 
for the Higgs mass m#, the critical temperature T c and the ratio <fi c /T c according to 

m H = 8Xv , — = — , T c = — = . (59) 

H Tc 2A' v/8 Cl A - 2c| 1 ^ 

The latent heat £ c is given by 

l c ~Td T [V(T, <j) c ) - V(T, 0)] T=Tc = 4A0;V , (60) 

and hence 

ac = ■ (61) 

The nucleation temperature can be determined using the expression for the tunneling 
action in the thin-wall approximation 

S3/T = V2 mHfV[ 1 -T c ) ' (62) 
yielding for the nucleation temperature (S 3 /T N w 140) 

log— 0.04— — \ — . (63) 



Using this in (14011 one obtains 




ee( 07,,4)(^)l-H |?W£, (64) 
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Figure 3: The wall velocity of the obstruction versus (f> c /T c for the values T c = 0.6, T c = 0.9, and 
T c = 1.5. The Higgs mass is m# = 120 GeV and the gray bands correspond to the estimate in ([40]) . 



where the number of degrees of freedom in the symmetric phase g* is given by g* = 30 a + /-K 2 . 
Hence the obstruction velocity depends in the regime of weak phase transitions only weakly 
on <f) c . Still, the thin-wall approximation ( 1621) is not very accurate and we use numerical 
results for the nucleation temperature in the following. 

The result from solving the local equations compared to the estimate in (j4"0l) is shown in 
Fig-El As expected, for small velocities the estimate works rather well while it underestimates 
the heating effect in front of the wall when the wall velocity approaches the speed of sound. 

In the last section we found the criterion ( I48p for the occurrence of the obstruction under 
the assumption that the temperature in front of the wall equals the average one. In the 
following we discuss how this criterion is modified if this assumption is replaced by the local 
system of equations. Still, one would like to express the criterion for the occurrence of the 
obstruction in terms of critical and nucleation temperature only, since these quantities can 
be obtained by inspecting the effective potential without an hydrodynamic analysis of the 
deflagration mode. Especially for weak phase transitions one would expect that the criterion 
(jlSj) is only modified by a proportionality factor c(a c ) close to unity: 

log ^ <c (a c Ma c ) . (65) 

In order to test this hypothesis, we determine the wall velocity of the obstruction in 
the model (155]) . We keep a c fixed and vary T c in such a way that £ w = c s . The functions 
log lp- = Q(a + ), log y~ and f2(a c ) are plotted in the right panel of Fig. HI It turns out that 
the proportionality factor in (1651) is for small a c close to |. This factor depends in principle 
on the model but it is expected to be in the range [|, l] . The upper bound results from the 
fact that the average temperature in the wall cannot exceed T + . The lower bound is due to 
the fact that the average temperature cannot fall below T_ and that for small a c the heating 
effect is behind the wall half as strong as in front (see eqs. ( 131)]) and ( 155]) ). Still, there is a 
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certain uncertainty in this factor^. 



0.25 



0.15 



0.05 




O log(T + /T N ) 
□ log(T/T N ) 



Figure 4: The heating effect assuming a wall velocity of the speed of sound, £ w = c s . The plot 
displays the relation between the parameters of the phase transition at critical temperature in 
relation to the parameters in front of the phase boundary in the toy model (|58p . 



Interestingly, the criterion (1481) becomes exact in the limit a c — > 1/3. This can be 
explained by noting that in the limit a c — > 1/3, the velocity in front of the wall equals the 
bubble wall velocity, v + = 0, and so the temperature behind the wall T_ has to drop to zero. 
Thus, the combination u^ 2 v is constant in the wall and small. According to (I5T1) this implies 

dcj) 1 OtP 

dT ' 



(66) 



and therefore the Higgs vev changes quickly close to the broken phase (where v is small) 
and slowly close to the symmetric phase (where p is small). Hence in this limit one obtains 
T + — > T c and the criterion ( l4"8"j) becomes exact. However, this behavior could result from the 
assumption of local equilibrium and might not be reproduced if the full set of Boltzmann 
equations are solved instead. 



5 Friction revisited 

The analysis in the last few sections dealt with the hydrodynamic obstruction in the limit of 
vanishing friction what gives an upper limit on the wall velocity. In this section we indicate 
how both effects - heating and friction - can be taken into account. 

Our starting point is again the equation of motion of the Higgs field f|T2|) 



□0+^-/C(0) 







(67) 



2 For a model with an additional 4> e term in the Higgs potential the proportionality factor also turns out 
to be close to | similar to the toy model of weak phase transitions. 
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where J-{4>, T) denotes the free energy of the system and JC(4>) quantifies non-equilibrium 
effects that ultimately lead to friction and hinder the bubble expansion. Integration of the 
equation determines the acceleration rate of the bubble wall, which is of the form 

4 oc (Jl - F+) (T) - £ wV , (68) 

where J-± are the free energies in the symmetric and broken phase, and the coefficient rj 
parametrizes the friction. The temperature T should be identified with the average temper- 
ature in the wall, that can differ significantly from T/v due to the compression wave in front 
of the bubble wall in the deflagration mode. Taking this effect into account leads (for small 
wall velocities to an equation of the form 

£ w oc (JL - ?+) (T N ) - i 2 w K - £ w r] , (69) 

where k is (for small a c ) approximately given by 



T J dT 



T=T C 



C ( 1 - Y- ) 4 * 4a c 4 . (70) 



This equation is only valid for weak phase transitions (a c <C 1) and small wall velocities 
(C,w "C 1) but a generalization can be easily obtained by solving the hydrodynamic system 
along the lines of sec. [3] in other regimes. For models where the friction coefficient 77 is known 
(in particular for models with a particle content similar to the SM or MSSM), this approach 
can be used to obtain accurate wall velocities. 



6 Discussion 

We discussed the expansion velocity of nucleated bubbles in a first-order phase transition. 
For the subsonic deflagration mode, the expansion is not only hindered by friction (that 
results from deviations from equilibrium) but also due to a heating effect by the compression 
wave in front of the bubble wall. Even if friction is neglected, the heating effect by itself can 
lead to an upper limit on the expansion velocity of the bubbles. 

In summary, a criterion for the occurrence of this obstructions at subsonic wall velocities 
is given by 

log ^ <c (a c )n(a c ) , (71) 

where T c and T/v denote the critical and nucleation temperature of the phase transition, 
respectively. The parameter a c quantifies the strength of the phase transition and is given 
by « c — 4 where £ c denotes the latent heat and e c the energy density of the plasma in the 
symmetric phase, both evaluated at the critical temperature. A fit to the mo del- independent 
function fl(a c ) is given in eq. (jl7)l . Finally, the function c(a c ) is a model-dependent fudge 
factor that approaches unity for a c — > | and a constant in the range [|, l] in the limit 
a c — > 0. For the model of weak phase transitions given in fl58l) as well as the SM with an 
additional 6 term in the Higgs potential it is found to be close to | . 

Let us discuss this result in the context of different models. Consider the MSSM in the 
light stop scenario [2HI2S]- Using the characteristics of the phase transition 

C ~ T c ~ 95 GeV, m H ~ 120 GeV , (72) 
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in ( 16"Tj) and (16"3]) leads to 



T 

a c = 2.6 x 1(T 2 , log -f = 3.4 x 10" 3 , (73) 

J- n 

and an obstruction at the velocity £ w ~ 0.18. Hence, one can conclude even without the 
knowledge of the friction coefficient that the wall velocity in the MSSM is subsonic what is 
a necessary requirement for electroweak baryogenesis. An analysis of friction without the 
heating effect [TT] leads to a wall velocity £ w ~ (0.05 -r- 0.1) what shows that in this model 
friction is very relevant. 

Next consider the SM with the (unrealistic) parameters 

C ~ 60 GeV , T c ~ 120 GeV , m H ~ 90 GeV , (74) 

what corresponds to one parameter set in ref. [10] and yields 

a c = 2.3 x 10~ 3 , log -f = 1.1 x 10~ 3 , (75) 

J- n 

and hence ^ w ~ 0.35. In comparison, in [10] the heating effect as well as friction was taken 
into account and a wall velocity in the range £ w G [0.35, 0.40] was found. The result was not 
very sensitive towards changes in the interaction rates such that the heating in front of the 
phase boundary is the dominant effect hindering the wall expansion. This also fits well with 
the results in [27] that neglected the heating effect and found much larger wall velocities. 

Let us discuss some cases where the wall velocity is still unknown. Consider the nMSSM 
as discussed in [281129] . In this model, typical parameters for a phase transition that is on 
the strong end of the spectrum ard^l 

T c ~ 75 GeV , T n ~ 65 GeV , a c ~ 8.9 x 10" 2 , (76) 

what leads to a subsonic wall velocity according to ( 171]) and electroweak baryogenesis is 
feasible in this model. 

Finally, consider the SM enhanced by a large number of singlets from a hidden sector 
as discussed in [25J. For a singlet coupling to the Higgs of order £ ~ 1.0 and a Higgs mass 
mH = 125 GeV one finds for example the following characteristics of the phase transition^ 

T c ~ 100 GeV , T n ~ 93 GeV ,« c ~1.9x 10~ 2 , (77) 

and hence a wall velocity £ w ~ c s [using the criterion for large velocities (ITT]) ]. For smaller 
values of the scalar self-coupling ( and/or larger Higgs masses the phase transition becomes 
weaker such that for these parameters the wall velocity is definitely subsonic. On the other 
hand, under the assumption that the friction is similar to the SM, the wall velocity is 
supersonic for larger values of ( and/or smaller Higgs masses and can even enter the runaway 
regime according to the analysis in ref. |14j . 

3 Notice that in j28j[29] the latent heat normalized to the energy density at the nucleation temperature is 
denoted a. 

4 Also in |25) the latent heat normalized to the energy density at the nucleation temperature is denoted a. 
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